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Abstract 
Phenomenological analysis to the R-parity violating supersymmetry with a vector-like extra 
generation is performed in detail. It is found that, via the trilinear couplings, the correct neutrino 
spectrum can be obtained. The Higgs mass rises to 125 GeV by new up-type Yukawa couplings of 
vector-like quarks with no need of very heavy superpartners. Phenomena of new heavy fermions 


at LHC are predicted. 
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I. INTRODUCTION 


Recently_a standard model (SM) Higgs-like particle with a mass of 125 — 126 GeV was 
discovered |1]. In the paradigm of the weak scale supersymmetry (SUSY) which aims at 
the naturalness of the electro-weak scale, however, such a Higgs mass brings in tensions, 
especially the minimal SUSY SM (MSSM). Nonminimal and still natural scenarios of SUSY 
are thus motivated. One of them is the MSSM with a vector-like generation Hu It gives 


has a rich phenomenology . In the framework of SUSY, vector-like fermions can also 


the right Higgs mass naturally, is consistent with precision electroweak measurements, and 
be motivated by other theories beyond SM, such as SUSY extension with extra-dimensions 
or with composite states DÉI. So it is worth asking the question whether such a scenario also 
provides explanations to other problems such as neutrino masses. 

Neutrino oscillations are the undoubted new physics beyond the SM. Daya Bay |8| and 
RENO [|9| experiments recently discovered a relatively large 013 ~ 8.8? + 0.8°. Within the 
framework of SUSY, in the absence of R-parity conservation, neutrino masses and mixings 
can be generated from lepton number violating (LPV) couplings ld. This approach was 
extensively studied before ll. It is known that all the neutrino experimental results, includ- 
ing that of oscillation phenomena like the large atmospheric mixing angle 055, the hierarchy 
of oscillation frequencies Am2, « Am, and the smallness of 0,3, can be understood in 
three generation LPV MSSM. However, this needs some special requirements for relevant 
coupling constants and mass parameters. 

Combining both considerations above, we will work in the LPV MSSM with a vector-like 
extra generation |4]. While this model takes the vector-like slepton doublets as the two 
Higgs doublets needed for the electroweak symmetry breaking, the SM-like Higgs mass can 
be naturaly 125 GeV [5]. Extra trilinear LPV couplings between ordinary fermions and 
vector-like fermions provide a much larger parameter space to explain neutrino pheomena 
right. 

In this paper, phenomenological aspects of the model will be analyzed. In Sect. II, we 
make a brief review of the model. In Sect. III, neutrino masses are calculated. For the 
neutrino physics, noting the enlarged parameter space, we consider trilinear LPV couplings 
carefully. One-loop contribution to neutrino masses due to new trilinear LPV couplings 


is calculated, theoretical analysis are performed and numerical results are shown in detail. 


Besides, we analyze the SM-like Higgs mass and explicitly show that it can be increased to 
125 GeV by two new Yukawa couplings of the up-type Higgs with vector-like quarks in Sect. 
IV. The LHC phenomenology of the new fermions is analyzed in Sect. VI. The summary 


and discussions are given in the last section. 


II. A BRIEF REVIEW OF THE MODEL 


This modelld is SUSY and SM gauge invariant, and R-parity violation with baryon 
number conservation is assumed. For the matter content, in addition to the ordinary 3 
generations (3G), a vector-like generation is introduced in. Without R-parity conservation, 
this can be also thought as that there are 4+ 1 chiral generations, where '4' stands for four 
chiral generations with SM quantum numbers and ’1’ for another chiral generation with 
opposite quantum numbers. The 4 chiral generations with same quantum numbers mix. 
The ’1’ has Dirac masses with only one combination of the ’4’, thus, there are always SM 
required three massless chiral generations and one massive vector-like generation. 

In terms of mass eigenstates (before electroweak symmetry breaking), the massive slep- 
tons in the vector-like generation are taken as the two Higgs doublets. New particles beyond 
the MSSM are the following with quantum numbers under SU(3),xSU(2),¿xU(1)y, 

Ej(,1,2) , E&(,1, 2) , Qa(3,2, 3), Qr(8,2,-5) , 
U¿(3,1,-5), U£(3,1,5), D&(3,1,2) , D&(3,1, 3). 


The superpotential is conveniently written as 
W=W,+Wz, (1) 
where Wo and Wy stand for that with lepton number conservation and LPV, respectively, 


Wo = 4H, Ha + y EGEj + un QíQg + YUGUR + uP DUDE + yl ¿Li HAS + y Hu De 
+y44Q¿H,U; + yP L;HGE; + VO, Har + yPQ Haf: + y??QaHaD§ + y Q¿H,Uz 
+y; QU? + y" QLEU$ + y Oo HAUS, + y "QuH,Dj,, 


and 


Wy 2 Aur, iL; ¡Ex SS Mi Qi; Di. + AG, Lj Es + AGO D; 
NPQ DI EAPPQUEDS EES, (2) 


where Li, Qi, Ef, Df, Uf, i=1-3, are the first three generation SU(2)z doublet leptons, 
doublet quarks, singlet charged leptons, singlet down-type quarks and singlet up-type quarks, 
respectively. H. and Hq are the up-type and down-type Higgs. Note that the term Q y HD, 
in Wo was missed in Ref. HIT And in Wy interactions of purely singlets are omitted, which 
are irrelevant to our study. 

By assuming universality of the mass-squared terms, the alignment of the B terms the 


soft mass terms and the trilinear soft terms of all fermion's superpartners in the model are 


-L D MIL, + ME Ha + MPH + MEEÑES, + AO Qm + MÜA, 
EMPDADE, + May EZ Eg + Mog Qu Qu + MogUgUg + Mäe Dë D, OG 


T(BuH4H, + Bew ESEG + Div, On + BU" U;Ug + B'uP DD; + h.c.). 


Proper values of the new DH Hund HH terms are set to avoid unwanted color symmetry 
and purely U(1)y symmetry breaking, see Eq. (11, 12) in paper [4], therefore EWSB in our 
model is just the same as in MSSM. After EWSB, the specific fermion mass matrixes and 


sfermion mass-squared matrixes are given in Appendix A. 


III. NEUTRINO MASSES AND MIXINGS 


LPV results in nonvanishing neutrino masses. In this model, in addition to traditional R- 
parity violation in the MSSM, a lot more bilinear and trilinear LPV interactions are brought 
in through the vector-like generation. In this work, the trilinear R-parity violating interac- 
tions will be studied. To avoid complication due to too many LPV sources, sneutrino VEVs 
will not be considered. There are several reasons for this.First, we can phenomenologically 
assume the universality of the soft SUSY breaking mass terms at the weak scale, to avoid 
dangerously large flavor changing neutral currents (FCNCs), without considering any UV 
completion of the model. In that case, because of the alighnment in bilinear terms of the 
superpotential and that of soft terms, R-parity violating bilinear terms can be rotated away 
via field redefinition, and sneutrino vacuum expectation values (VEVs) vanish in the phys- 


ical basis. The second reason is from consideration of underlying models. SUSY breaking 


! It modifies the down-type fermion mass matrix and scalar mass-squared matrix. Correct ones, as well as 


the resulting mixing matrix are given in the Appendix A. 


is introduced effectively in our model, it can result from gauge mediated SUSY breaking. 
Then the messenger scale can be as low as 100 TeV, even if the universality scale is at the 
SUSY breaking messenger scale, the running effect is small, and the bilinear LPV is not 
important compared to the trilinear ones. Finally, small sneutrino VEVs can be included 
in the analysis nevertheless in future works, after the role of new trilinear LPV interactions 
gets a thorough understanding. 


The trilinear LPV Lagrangian relevant to neutrino masses is from Wy, 
L C age (linPinlin + lirlervir) — Nie (dkrUindye + dirdervir) 
Ni (Bibir + LEY vi) — Ay (digPinQa + Hadamar) (4) 
—AF (O50 Qn + QuUg vir) + h.c.. 


where 775 stands for the left-hand neutrino. 

The 7 types of trilinear LPV interactions in the above equation induce 14 types of one-loop 
diagrams contributing to the neutrino spectrum, which are proportional to AA, A’, APA, 
AXE, ARAM, ACAP, APAP, AMAM. ARPA, NAS, NAP, NASP, A AN respec- 
tively. The Feynman diagrams and the corresponding analytical results are shown in Fig. 1 
in Appendix B. For simplicity and without losing our purpose, in the Yukawa interactions of 
Wo we assume that only af. y9', y9', yP, y", y”, y" are nonvanishing, that is vector-like 
particles have Yukawa interactions only with the third generation. Thus, the vector-like 
generation has little constraints from the collider phenomenology. 

Before starting to analyze the neutrino mass spectrum, some assumptions are intro- 
duced in order to control the parameter space and get relatively simple analytical result. 
q 


Since four new up-type Higgs Yukawa couplings y", y9, y9", y" and five new down- 


type Higgs Yukawa couplings yP, yP, y9', y9P, y” 


appear in our model, and among 
which y9P, y9U, yĦ, yf provide the mass mixings between vector-like generations, and 
further more, they have infrared quasi-fixed point " so we assume y9P = y = 0 and 
y" ~ y = yt, < 1. We also set y? = y% = 0, y? < 0.04 and y" ~ y? = yl, < 0.08. In 
other words, we neglect all new down-type Higgs Yukawa couplings in quark sectors while 
consider all of the new up-type Higgs Yukawa couplings only and take yi, < y\,, which is a 


reasonable assumption. 


Basing on the above assumptions, contributions from AA, AY type diagrams can be 


simplified to the familiar forms 


TR 
y : 
ME S — 5 Ai33Aj33 Mr Sin Az cos oz; In —, 
8T TL 
M| > (Aisa; ino; sl br Aa, in a; zl br 
ij XX = gm Vas ¿334 SIN a, COS Qr IN — + Aj234532Ms SIN AG COS az In — 
T br br 
SR 
+ Aus 23 Mp sin a cos o; In —), (5) 
SL 


where in the first equation, we only keep the dominant contributions and in the second 
equation, we keep the dominant and subdominant ones. a@,, oy, Œs, o; are the angles of 
the corresponding 2x2 Tra», drin), Sum, fum unitary matrices. Unfortunately, the other 
equations, (B3)-(B14) in Appendix B, can not be simplified by following similar process 
because there are mixings between different vector-like generations. So these can only be 
analyzed numerically and will be discussed later. 

At last, without loss of generality, among all 7 types of LPV trilinear couplings, we take 
4 of them, AE, AP, AN and A”, for consideration while assuming the rest of them, A, A and 
AND. are negligible. The realization through different LPV trilinear coupling combinations 
can be derived straightforwardly. The method to calculate the neutrino mass matrix we use 
is given in Appendix C. 

Here we list the parameters of neutrino oscillation given by experiments, Am}, = (7.59 + 
0.21) x 10-5 eV”, Am?, = (2.43 + 0.13) x 1073 eV? and sin? 2045 = 0.861* 0025, sin? 2093 > 
0.92, sin? 2013 = 0.088 + 0.008. Scanning the parameter space with proper EWSB, we find 
by adjusting the ratios and values of the LPV trilinear couplings we choosing, the correct 
neutrino spectrum can be generated through the AFAF, APA9 and AMA? type one-loop 
diagrams. Numerical illustration is shown in Table I, in Set I we take the mass mixings 
assumptions mentioned before, in Set II we take different mass mixings and bigger vector- 
like masses for comparison. The specific parameters settings see Appendix C. 


That is by choosing (for Set I) 


AT AS 
0, 214, 5921x105, Ag M, ABW~ABs A A 
A23 A33 
A Aj H Q E E E Q 
rä, uml. An Ada Ata An Y Ags A, (6) 
2 2 
we have 
Mya —4 —6 
e 0.17, Mma ~ 5.9x 10 ^ eV , m,ı »5.1x 10^ eV, 
Mp3 
sin 013 ~ 0.143 ; sin 055 ~ 0.581 ; sin 01 ~ 0.559. (7) 


0.0027 0.168 0.004 0.003 0.011 


TABLE I: Numerical illustration for 5 types of one-loop contributions in our model ,the specific 


parameter settings see Appendix B. M7; (GeV) stands for the parts in Eq. (4,5) excepting the 


LPV trilinear coupling constants. 


Unlike in the 3G LPV case, where Ajz3A533, Aj934532 + AjgoXjo3 and Aa Aa type one-loop 
contributions are dominant, subdominant and next-to-subdominant, here in our model, un- 
der the assumptions mentioned before, Ee $55, AR Ad and AA type one-loop contributions 
are dominant, subdominant and next-to-subdominant, respectively. This is because the new 


fermions 71, Dua, b12 in the internal lines, see Fig. 1, are much heavier than the third 


generation fermions 7, t, b. 


T,T1 b, De t,t12 
PEL unes ET 9 
E / NAE D / \ \Q H / \ AH 
Aj Z j A33 Ac " V AS A; 1 i Aj 
> L + + > L + + > L + + 
Vi V; Vi V; Vi V; 
iL PT jL iL b by B jL iL t, t12 jL 


FIG. 1: New one-loop contributions to the the neutrino masses and mixings from AFAF, AQAP 


and AP A type couplings. All particles stay in mass eigenstates. 


For the same reason, our requirements of the new LPV couplings we choose are of order 
10~° and small enough to avoid measurable FCNC decays such as u — ey . It worth 
to note that by decoupling the vector-like generation, correct neutrino masses and mixings 
cannot be obtained via AA, AA type one-loop contributions. 

In addition, AN AN type contribution containing up-type (s)quarks in the internal lines is 
absent in 3G LPV models because the vector-like down-type doublet quark Q^; mixe with 
the right hand singlet top quark. 

From Table I, we can also see that by choosing e SA, AN ye and AB Aja type one-loop 


contributions, the correct neutrino spectrum can also been generated in parameters Set II, 


we don't list the detailed results here. 


IV. HIGGS MASS 


There are four new up-type Higgs Yukawa couplings in our model, y", y9, y, ul, 
corresponding to the Yukawa mass, m§,, miz, m4, m',, separately, and also five new 
down-type Higgs Yukawa couplings, af. yP, y9', y9P, y”, corresponding to the Yukawa 
mass, m34, M34, Mz, m, mij, separately. The related superpotential contributing to the 
lightest scalar Higgs mass is shown in Wo. According to the assumptions mentioned in the 
last section, we neglect the down-type Higgs Yukawa contributions and the small up-type 
contributions between the SM third generations and the extra vector-like generations. The 


relevant superpotential can be simplified as 
W = 12Q1Qu + D Die + y" Qu D, D + yQ HU, , (8) 


So when neglecting the small D-term and the two-loop contribution, the new one-loop 


contribution to the lightest scalar Higgs square-mass is |5, |16 


3x2 1 1 1 X? 1 
Am = > Wt sinte ftv — ¿(5 0-a O 


where v = 174Gev indicates the Higgs VEV and 
H' U 2/472 M3 
y? =y? =yy , x = Ms/My , tv = logy (10) 
V 
(Au — Ham cot 8)? = (Agu — Lau cot 8)? = (Ay — py cot B)? = Xy, 


in which, for simplicity, u9 = uP = My stands for the vector-like mass of the new up-type 
quarks, M = Mj, = m? (see Eq.(4)) and Ms = vi Mg + m? stands for average mass of the 
new up-type squarks. 

In MSSM, the Higgs mass from the t, t one-loop contributions is about 110 GeV, for 
At = u = 400 GeV, m; = 400 GeV and tan = 10. Direct search bounds from CMS 
for exotic heavy top-like quark set limits of My > 557 GeV if B(t' > Wb) = 1 and 
My > 475 GeV if PI — Zt) = 1 [18]. When considering the mass mixing between the 
vector-like quarks and the SM third generation quarks, in other words, considering the 
realistic branch ratios, the mass limit is adjust to be My > 415 GeV bd. 5o if we set 
the vector-like fermion masses in our model to be My ~ 500 GeV, the soft supersymmetry- 
breaking parameters to be m ~ 700 GeV, Ay = uy ~ 500GeV and By uy ~ 500? GeV?, then 
from Eq. (10), in order to get approximately 125 GeV Higgs mass, for about My — 500 GeV 
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and Ms = 850 GeV, we just need to set yf, ~ 1, or say, need to set mi, = mo, = mí, ~ 174 
GeV. These values are just near their infrared quasi-fix point, as mentioned in last section. 
Evoked by the ATLAS and CMS discovery of the enhancement in yy channel and little 
deviation in ZZ channel Bi A the effects of the exotic vector-like quarks to the Higgs 
production and decay have been extensively studied recently . In general, in a theory 
with N vector-like generations extension, the new fermion contributions are suppressed by 
N? mi? / M? " . So only the very large couplings to the Higgs can obviously enhance 
the Higgs production and decay in the yy channel D but as we have mentioned, these 
couplings have quasi-fix point which limits their TeV values to be about 1 [5]. This value is 
large enough to accommodate m, ~ 125 GeV, but too small to influence the Higgs decay, 
one can’t depend on vector-like fermions by themselves to modify the Higgs decay branching 
ratios. As far as the Higgs problem to be concerned, extra vector-like fermions are mainly 
introduced to adjust the Higgs mass. However, the yy and ZZ channel anomaly, if they 
"i which beyond our scope in this 


persist, can be realized through the light stop scenario 


paper. 


V. THE EXTRA VECTOR-LIKE FERMION DECAYS 


To be clear, we list the new extra vector-like fermions below: 


c t,b c c 
= e Bez ut o Uy És Ms i , (11) 
E E Ug Dg 
in which E% mixes with zu: Ej mixes with rg; Q}, Df, mixes with bz; D$, Qi; mixes with 
br, Q4; Uf, mixes with tz, U$, Q^; mixes with tg. These exotic heavy fermions can decay 
into SM bosons, see Fig. 2, which will analyze bellow. Our analysis agree with the results 
given in [5]. However the slightly difference comes from their neglect of the contributions 
proportional to sz, in the vertex of Feynman rules. 
Note that theoretically speaking, when kinematically allowed, the exotic fermions pre- 
dicted in our model have the other two decay modes: through supersymmetric gauge kinetic 


interactions or the supersymmetric Yukawa interactions, decay into chargino/neutralino and 


sfermions, such as rı > C+v,, by > Njb, tı > C-b , where N;, i=1-4, is neutralino, C* is 


chargino; through LPV interactions, see Eq. (2), decay into fermions and sfermions, such as 


T1 T T1 i ` T T1 g ` U 

jp Pd Ka 
bio b, by bio b, by br tt 
the tty tie bti tie b, by 


FIG. 2: Tree-level decay of new exotic fermions in our model, all fermions stay in mass eigenstates. 


Tı > eji, bı > Tt, tı — 6b. Although the kinematical conditions for the latter decay mode 
are easy to be satisfied, but we have already seen in section III, the LPV couplings in our 
model, in order to explain the neutrino spectrum, are of order 1079, so we can neglect this 
kind of decay channels reasonably. On the other hand, for simplicity here in our work, we 
assume the former decay mode is not kinematically allowed. Therefore, the exotic fermions 


can only decay into SM bosons. 


A. T7; decays 


The weak bosons interaction Lagrangian to T, 7, is 


W = m — Z = m Z = m 
LO Ihr TIL VLW, + In EY TL + Uns TIRY TRZ p 


o L @ 
95, rg Del T S nm + h.c., (12) 


the couplings and the decay widths of 7; are given in Appendix D. 

The main characteristic of the lepton sector is that there must be mass mixing between 
the third and the vector-like lepton, otherwise the new heavy charged leptons 7, will be 
stable and give unacceptable cosmological heavy relic bd. For Specific, when y* — 0, the 
off-diagonal elements of L7, RT are equal to zero. That's why we set y" Æ 0 in section 
II while discussing neutrino spectrum, more specifically, we set y" < 0.04. Under these 


parameters settings, numerical results of 7, decay into W, Z,h? are shown in Fig. 3, we 
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FIG. 3: The decay widths of the new lepton 7; (left panel) and its branching ratios (right panel) 
with y? = 0.04. 


can see in the limit of large m,,, the branching rations are Bäim > Wv,) ~ 0.7 and 


BR(r; > Zr) = BR(n > hr) ~ 0.15. 


B. 11» decays 


The weak bosons interaction Lagrangian to t, ti, ta is 


L3 gi. a, try b, W,, 4 9.0, Por b, W + gi. os luv bRrW,, + gr. as DRY bel 
Ir Pr PW, + ges ton" iW, + Sun hi" tu Bu + Sun bv" tZ, 
+9 n oor" Zp zb IE as hi Y ERZ B Jé s Rl EZ F JÉ ntin ER" tırZu 
Jo p r 1 af, a ra aa 


+91 „ortırh” + ge, hirta? + h.c., (13) 


the couplings and the decay widths of t; ? are given in Appendix D. 


As mentioned in section IL, we take y" ~ yl < 0.08, y" ~ ul < 1, the numerical 
reasults are shown in Fig 4, 5. We can see in the limit of large Mi, ,, the branching rations 
of tı are BR(t; > Wb) ~ 0.4 and BR(tı — Zt) = BR(t, >H%) ~ 0.3, the branching 
rations of t9 are BR(ta > Wb) ~ 0.85 and BR(ta > Zt1) ~ 0.15. 


11 


0.020 r r r r r r x 


0.8 + 1 
3 
0.015 | b: 
\ 
$ E 
E 0.6 | \ 
= x 
= L E s. 
= 0.010 E So 
a 
Ki 
Q 
v 
a 
= 


0.005 + 


0.000 0.0 
300 — 400 500 600 700 800 900 1000 300 — 400 — 500 600 700 800 900 1000 


Mu (GeV) Mu (GeV) 


FIG. 4: The decay widths of the lightest new up-type quark t (left panel) and its branching ratios 


(right panel) with y@? = ul = yP — 0, yU ~ y = 0.08 and y" ~ yt — 1. 
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FIG. 5: The decay widths of the heaviest new up-type quark t» (left panel) and its branching ratios 


(right panel) with y@? = ul = yP =0, y" ~ y = 0.08 and 4" ~ y = 0.98. 
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C. bj» decays 
The weak bosons interaction Lagrangian to b, bi, ba is 
LD gp, LY tW, + of, bory tW, + gg pory ERW,, + gp, us Don" ERW,, 
Lo 1, LV LW, + gg uu DRY RW, + gg s bur "DL Zu + gj, pr Do Y OL Z 
+96», Don Y iL Zp, T Te, o DRY ORZ y + o p, LR Y ORZ + JÉ a, LR Y dırZu 
Ü = 0 yu o 7 o = 
+95 i Publ? + 9,1, ,bcbinh® + gr, Bahn + 9È p br boh? 


+ un barburh" + DEEN + h.c. (14) 


the couplings and the decay widths of bı 2 are given in Appendix D. 
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FIG. 6: The decay widths of the heaviest new down-type quark bg (left panel) and its branching 


ratios (right panel) with y9P = y = yP — 0, y" ~ 9 = 0.08 and y9" ~ all = 0.98. 


The numerical results of bə decay widths and branching rations are shown in Fig. 6 
under the parameter settings mentioned before, we can see BR(bz — Wb) ~ 0.85 and 
BR(b > Zbi) ~ 0.15. The branching rations of b; are BR(b; > Wt) = 1 which are not 


shown here. 
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VI. SUMMARY AND DISCUSSION 


We have studied several phenomenological aspects of the LPV MSSM model with a vector- 
like extra generation: the neutrino spectrum, the Higgs mass and the LHC phenomenology 


of the new predicted fermions. The results are: 


e The correct neutrino masses and mixings, especially the relatively large 0,3 can be 
generated from trilinear LPV couplings. The new trilinear R-parity violating couplings 
make it easy to generate the proper value of 0,3. These coupling constants need to be 


about 107*. 


e The two new up-type Higgs Yukawa couplings, al" and y9", between the vector-like 
quarks and the SM third generation quarks, with values about 1 near to their infrared 
quasi-fixed point in TeV scale, can give rise to 125 GeV Higgs mass with no need of 


very heavy new superpartner. 


e There are five new heavy fermions, 71, t1», b1 2 , predicted in this model. They can only 
decay into SM bosons by some kinematic assumptions. The branching radio depend 
on the mass mixing between the vector-like fermions and the SM third generation 
fermions. These charged exotic fermions would be quasi-stable if such mass mixings 


are very small. 


Based on our previous work about bilinear LPV couplings, further research on the renor- 
malization group (RG) of them is worthy to be studied in the future. There are also plenty 
of aspects to be further analyzed in the area of new fermion LHC phenomenology based on 


this model. 
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Appendix A: THE (S)FERMION MASS MIXING MATRIXES 


Because the CP violation is not considered in this paper, we have taken all the masses 
real. In this model, the mass matrix M of the third generation lepton and the vector-like 


lepton is given as following 


Lo -(, EM | |, (A1) 
Ez 
and 
Mt = M33 M34 (A2) 
0 u? 


U v — 
where m$, = Ys 7 cos B and m3, = y —= cos B. Taking |u®| > |m3,|, then the biunitary 


v2 


matrix to diagonalize M7 are 


L™ M Rİ = (mr, u”) = diag(m,, Mri) : (A3) 
where 
"Rx qw 2 (A4) 
u Tb LL E 
H (uE)? 


The mass matrix M? of the third generation down-quark and vector-like down-type 


quarks is given as following 


pe 
Lo- (b, Q}, Dg) M? | Dei, (A5) 
Qu 
where 
más má, 0 
M’= | mb, mb, u? |, (A6) 


0 pP mh 


v 


v2 


Taking that |u®| ~ |u?| > |m3,|; [mia], Imbal, [m3,|, then the biunitary matrix diagonalize 


M! are 


v U v 
b — d b — yD b = ya b OD 
where m34 = ie cos B, m$4 = ya 3 cos B, mis = ys 5 cos B and mj, = Y cos D. 


L* MOR = (mp, u9, u?) = diag(m,, mg, mi) , (A7) 
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where 


1 0 2 
Dmt tum? 
P=) 0 1 oa | (A8) 
me (m? (m5 y 
Pa ane Fu m I 
and 
1 ma 0 
b uo matu? m} 
R= | 0 Pra? (A3) 
_ mba 1 (mir)? (m3)? 
un pP mi tu mj, 


The mass matrix Mt of the top quark and vector-like up-type generations is given as 


following 
ge 
£ 5 — (t, Qi, Uf) M' | UE 
H 
where 
M33 m$, 0 
A = mis má, UÊ | 
O uu me 
where mA, = o sin B, m$, = Wo sin B, mis = WF 


m^, = cos 3. Taking that |u9| ~ |un"| > [miz], mi, 


v2 


biunitary matrix diagonalize M‘ are 


(A10) 


(A11) 


U 
sin 8, mí, = y9" — sin B and 
44 — y Ka 


|, [más], [mái], [m |, then the 


LA RH = (m,, u, Y) = diag(m,, ma, ma) , (A12) 
where 
E 
1 0 — er 
t — uU ml tue ml 
dis RK TI | > 
mí, (mig) (má? 1 
uU mp mi, 
and 
1 ra 0 
t uQ mi, 4 uU mt 
SE (Gy Fu - (n : | ee 
E t 42 t 32 
EC à (máy (mb) 
DN p^ mie Eu mi, 


The charged slepton mass-squared matrix M? of 7 and the superpartners of the vector- 


like leptons is given as following 


Ls 
£35 Ey Es Es Eg) xe | I |, (A15) 
Ei 
Ej 
where 
(Min = M? + KE = HA cos28 + m2 + [mia]? , (MP)12 = (mo — u tan B)m,,, 
(M2)ıs = (mo — u tan B)m3, , (M3) = um , (Ma = (mo — p tan B)m, , 
(M?)o0 = Mg — (mz — miy) cos 28 + m? , (M?)o3 = momia, 
(M?)oa — 0, (M2?)31 = (mo — u tan B)m$, , (M2)32 = momia, 
(M?)33 = |" |? + ME — (m3 — miy) cos 28 + [mia]? , Lia = Bu, 
(M?)a =u myu, (Ma = 0, (Mas = BE uF, 
(A2) = |u?]? + M24 + (mZ — m2,) cos 28 . 
(A16) 
The corresponding unitary scalar matrix is defined as 
V* MV" = diag(M?, M?,, MÀ, MÀ), (A17) 


The mass-squared matrix M2? of b and the superpartners of the down-type vector-like 


fermions is given as following 


£> (8, Ds, Di Dg ÂY, 04) ME] |; (A18) 
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u =M¿- mari cos 28 + mi + m5? , (MP) 12 = (mo — u tan B)mp, 

13 = (mo — utan 8)m3, , (Mira =P mb, (M})ıs = memi + MEM 

16 =0, (M?)o1 = (mo — utan B)m, , (M?)o9 = M2, — má" cos 28 + m? + bag, 
23 = mom, + Mega, (Mijas = 0, (Mas = (mo — utan B)m2z , (Mi)as = miz, 
ai = (mo — utan lan, , (Mas = mum}, + migmiy , 

s = Pl? + MB — TEZE cos 28 + [mba]? + Imt? , (Mos = -BP pP, 

(mo — utan B)m34 , (Mi)as = mb, + Pm, (Mar = Pm, , (Mia = 0, 


2. m2 
mz my 


a = —BPy? , (Ma = |e? |? + Mën + TES cos 28 + mp? , 


a5 = Mim, + win, (Mi)as = (mo — uicot B)m', , (M3)51 = mimis + migas, 


52 = (mo — utan lan, , (Mi)ss = (mo — utan B)mhy , (Miss = ui? ml, + im, 


m2, m2. 
ss = [ul + Mô ES DTP. cos 28 + mi? + mis)? , (Ms = BERS , 
eu —0, (M?) 62 = um). > (LZ le E umi, + uP m, > (Moa = (mo — cot ea 


= m2, +2m? 
e = BI, (Mies = uel? + Mën + + mb? + FETTE cos 28. 


ee 
SN TN SN SN SN SN SN SN SN SN SN ON an 
—— AA A NA NS — u ut I A NL SA NS 
w 
or 


TM EOS OS us emm emm ES emm emm emm eme emm em 


(A19) 


The corresponding unitary scalar matrix is defined as 
VÄAGVN — diag( M7, Mà, Mb, Mé Mé Mi) H (A20) 


The mass-squared matrix M? of f and the superpartners of the up-type vector-like 


fermions is given as following 


ta CaaS Ue, Os, 020%) M? K , (A21) 
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(MP = M3 + T cos 28 + me + [mS P, (MP)12 = (mo — p cot )m,, 
(Mi) is = (mo — ui cot 8)m$, , (M2)i4 = H mi, (M2)is = mimi + MEM, 
(Mijas —0, (M2)n = (mo — u cot B)mi , (M?)20 = Mg + ¿(m3 — miy) cos 28 + m? + |mig|?, 
(M?)o3 = mimi, + migmiy , (Mia = 0 , (M2)as = (mo — pi cot B)m£s , 
(M2)as = miz , (Mia = (mo — u cot B)mb, , (M2)sa = mimi, + migmiy, 
(M?)33 = |u" |? + Mà + ¿(m3 — miy) cos 28 + [mia]? + mial? , (Misa = —B"u", 
(M?)35 = (mo — u cot B)mi, , (Mao = um, + umi, , (Ma =p msg, 
(MPa =0, (MPa = —BU pl, (M2) = |uU |? + An — 2(mZ — mi) cos 28 + |ms |? , 
(Mas = uU*má4 + uml , (Mas = (mo — utan B)m'y , (M?)51 = mamás + mama, 
(M?)52 = (mo — cot B)miz , (M?)s3 = (mo — cot B)mig , (M2)s. = p má + Hinne, 
(Miss = [u9]? + M3 + Bei cos 28 + miu? + misl? , (Mio = — Pin. 
(MPa = , (Me = Omis , (M)ss = u€ má, + uU mf , (M2). = (mo — u tan B)mi, , 
(Mies = —B9u? , (MP og = lf + Mg + Imi? — RP 00828. 
(A22) 

The corresponding unitary scalar matrix is defined as 

V¿M2V* = diag(M2, M2, Mà, MÀ, Mà, M2), (A23) 


Appendix B: NEUTRINO MASSES IN OUR MODEL 


All fourteen types of one loop Fyenman diagrams which can contribute to the neutrino 


mass and mixing in our model are shown in Fig. 1 
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FIG. 7: One-loop contributions to the neutrino masses and mixings in our model. 
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The corresponding analytical results are listed below: 


Y 1 *T DRT TYXTT/T 
Mila = [9 Kë Ai33Aj33 Rmi Lira Vet Vg Meri Das, Mz, ay): (B1) 
k,m 
V 1 XT *T *T T 
Mil = 82 A Aug Lo MET Vk3 ma, OM > Ms y) 
k,m 
+ RD Ver Vi Mr b(m;, Melo (B2) 
1 *T *T *T Y" 
Miles = 87? > Aja a Leni Vki Ves Mer OM, > Moos (B3) 
k,m 
v 3 * * * 
Mi be = 872 Y, Aus Aen Lina e my, DM y, Ms) 
k,m 
+ Aug Aan RO LO m. sin Qs1(2) cos a1 (2)b(ms,, , Ms, ll. (B4) 
Li 3 * * * 
Mijxexe = sr > ABAS Fo Lomo Vez Vis mn, HU Mis (B5) 
k,m 
v 3 * * * 
Mi ADAD = re ASA. Hr E Ve Vie Mom (Mp, Mi, s) (B6) 
k,m 
v 3 * * * 
Milano c 55 3. ANTANT Bri Lima Vis Ves mos (mo, Mi 7,1) (B7) 
k,m 
v 3 x * * 
Mr AQAD = Sa? p» ASA [Ro Lia Vei Vio Mm O(Mbm 5 Mean) 
k,m 
T Br Ua Ves Ves Mom b (mo, > Mi, (s) (B8) 
v 3 * * * 
k,m 
S Brno Dino Veo Ves Mom OM mn» Ai, (B9) 
v 3 * * * 
Milon E 872 >> AU [Ro La Vis Ves Mom (mi, , My os) 
k,m 
+ Bro La Vir Via Mom D, My isl (B10) 
v 3 f * * * 
Milne > es? b» AA Ri Dina Vez Ves m (Mb Mi, cae) 
k,m 
+ ER Mom (ms, , Mi ish (B11) 
v 3 ` D ab pb 37xb37b 
Mily 1D = En? > EA IN Vir ME Mm DA Ml 
k,m 
i Hr Usa Vin Ves Mm (mo, , Mi ish (B12) 
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v 3 L D * * * 
Milon = gc A. Ai} [Pro La Vez Ves my, OM Marima) 
k,m 
oF Rr Lra Ver Ves Mom b(mo,, , Mi, cay dl (B13) 
3 


Milan ~ bonos (B14) 


H\H pst rt *uY7zu d 
2 AT AF Rms Lms Vio Via me lt, M; 
k,m 


81? 


In which L^5', R ®* are biunitary matrices of mass matrices between (r,b,t) and the 
vector-like fermions (see Appendix A), while m.,,m»,,,m.,, indicate the corresponding 
mass eigenvalues. vit are the square mass mixing unitary matrices of their superpat- 
M5, 


ners, while M; 


BN Mi, Sat stand for the corresponding mass eigenvalues. sin Qs1(2), 


(ok? 
COS Qs1(2) are the unitary matrix elements of 5. b(mı, m2) is the loop integral factor: 
b(m,,ma) = ——— (m? In m? — mill ml — m3 + mi). The value range of the indices in 


Mima 


Eq. (4)-(6) is m=1,2, k=1-4 ‚while in Eq. (7)-(17), it is m=1,2,3, k=1-6. 


Appendix C: NEUTRINO SPECTRUM-CALCULATING METHOD AND PA- 
RAMETER SETTINGS 


The methods to generate neutrino masses and mixing angles with one-loop trilinear L 


couplings actually involves the following three matrices 


a? ab ac d? de df g? gh gl 
mı | ab b? bc | , m| de e? ef | , ma| gh h hl | , (C1) 
ac bc c df ef f? gl hl P? 


where we name each of the matrices above M123 separately We assume mı > mas, 
m» ~ ma and there is no strong hierarchy between a, b, c, d, e, f, g, h, 1. 


Mı has only one eigenvalue after digonalized by an unitary rotation 
XTMıX = diag(0, 0, Mi) , (C2) 
where 
M; = mi(aà? +P +e), (C3) 
and 


Ca S2C3 8283 
X = 523 C2C3 S283 > (C4) 


0 — $3 C3 
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a C 


AM === 
Var"? Varı+b+c 


$9 = 


If we rotate the sum over M1 23 by matrix X, it becomes 


€11 €12 €13 
XT (Ms + Mə + M3)X ~ mila +b? +e) | e e2 ex | > 


€s1 €32 1 


(C5) 


where €;; are some small values related with ma/mı, ma/mq and the other elements of Mj 23. 


We can then define another unitary matrix X' to diagonalize the matrix in Eq. (B6) in an 


approximate way: 


where 


and 


A CAM + Mə + M) XX’ RI ma (a? +b + e) diag(ö3, 05, 1) > 


Cy 51 0 
X = —9$1 Cj 0 > 
0 01 
2 
tan 20, = M . 
€22 — €11 


Then from Eq. (B7), we get all three mass eigenvalues 


Mı ~ m (a? +b? +e) ; M» ey M5 0, ; M3 os Mô; ; 


and from Eq. (B4, B8), we get all three mixing angles 


— — a 

$18 — $283 = VERE 
— b 

$23 — €283/C13 — Veg 


$12 = (5102 + C18203)/C13 . 


The parameter settings we used in table I are given as following 


Set I: 


mi, = 10, Me = Men = 600GeV, B*y* = 400°GeV?; my = 170GeV, 
má, = mis = Mig = 0, Mg = Mp = Mpx = Mon = 700GeV, 


B?P y? = Du = 500°GeV?; mi, = mi, = mi, = 13GeV, mi, = 174GeV, 


My = Mun = 700GeV, BY u” = 500?GeV?, tan 8 = 10, A = u = 500GeV. 
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(C7) 


(C11) 


Set II: 


m3, = 10, Mg = Mey = 1000GeV, B® u” = 600?GeV?; mà, = 170GeV, 


mi, = ml, = ml, = 10GeV, Mg = Mp = Mpy = Mon = 1000GeV, 


BP y? = B? ue = 600°GeV?; mi, = mi, = mi, = 13GeV, mí, = 174GeV, 


My = Mug = 1000GeV, Du! = 600°GeV?, tan 8 = 10, A = u = 600GeV. 


Appendix D: EXOTIC QUARK AND LEPTON COUPLINGS TO W.Z,h’ AND 
DECAY WIDTHS 


The couplings for the W, Z, h? boson with leptons in Eq. (12) are 


aus = ih , 
E = Pi ELI, = ios = Asiy)L5,L7,| , 
fuso SE RGR + B si RR) (D1) 
B us = E + y3 L5, Rip) > 
ho Sa 


= (y33 L4, Ra + y S Rao) , 


Irura EU 
where c, = 58, Sa = —cg is the elements of the rotation matrix related with the real parts 
of (H9, HW). 
Then the decay widths of 7, are 


Din > Wor) = 22.0 + ay — 22) PC a + tw) P. 
Dueb Sla =e, 95 9:923) 

E EE a M) T a 

+ 122,95, Dësen) > (D2) 
Dia > h?r) = TAN + rfo + xi — 222, — 222 — 232,52)? 


0 0 
(o 22 — lr)? + (ina)? 
h? he 
T 42:92 rari rn] H 


where x; = m;/m,, for i = W, Z,7, h?. 
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The couplings for the W, Z, h? boson with t, tı, t2 in Eq. (13) are 


W g b b W 
Ütpb = ell LU T 155161) ^ Üügbg — =R; itis > 


v2 3 
g 
Bä = vg Ena Y LG ’ Donde = SZ Sala E 


Ww g g b 
Üt,bn "7 N LE ES L5) > GE = vg Fan Ras > 


d sto + O — q (LL + Hoa) > 
Irin ch (Pia Ria + Ripa) + (2— sel) RR] 
ds, = E Is + Lo — SU + D (n3) 
ds = Te Get (RRi + RIS) + (2— Fsh) Ris] 
d = TRE IS + LO slo Us Us + Hale) 
GE ntin = E [$s (2$ ¡Roy + RS RS) + (2 — = Sty) Pis] ; 
P us — uis D Ri + ys La Ri, + Y Lo Fo + ya Lig Rig) — a y” Lais. 


Be 
Sa 
4/2 
Sa 
Wold 
Sa Hrt pt 
—y Lis R33 , 
ve 134233 
BY 


So 


v2 


GE Ltr v2 


0 Ca u U 
hr = sl us Un in 3r ys Dis + ys qmm gs y? Li; R$) 


de 


0 Ca QU 
JÈ tr = < (y5 L4 Ri, + y Leaks + ys ES + y? Lih aj — 


Ro 


0 Ca u U 
9E ton —Ñ(y33 Li R3ı KS URS + Us S Du Rio Es ys Lia Rio) z 


v2 


0 Ca u QU 
Meta = — (3305, Ry + UTR + Y3 dr + ys LR) 2)— 


V2 


0 Ca U 
GER = —(y33L21R3 jp Qu. + us Dus T y? L R52) — 


y2 


The decay widths of the lightest new up-type quark t; are 


yD igh , 


t 
eee Ris D 


y” t 
153R53 D 


y^ L5 : 


T(t; > Wb) = m "du qu O oa 
(+a? — 208, + afetten 
E Ir} > 

T(tı > Zt) = SCD +2 + r4 — Oe — Eier 


Ire — 222 + (1-27) z^) us + Lef en (D4) 


Z Z 
+ 122,95, 4, Inte) > 
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T(t, + ht) = S (1 + ao + of — 222, — 222 — 272,42)? 
T 
0 0 
{+ — ol + (Shan) | 


0 0 
E ? 
where z; = mj/m,, for i = W,Z,t,h?. The heaviest new up-type quark to has six decay 


channels. The decay widths have the similar forms and can be deduced straightforwardly. 


The couplings for the W, Z, h? boson with b, b1, b2 in Eq. (14) are 


g b b W 9 pb 
Put = ~za Li xd 13275) > onis = —Ro3Ri3 ; 


v2 a 


OE ty = ND T L5314,) , OE tn = KC > 

a due atb, + IL) > Intin = et SCH 

Ind = L33L43 — zl i Sate, + HEH , 

d = qe geh (Rl, + RIS) + (2 GR 

D d = P DIA, - do IC = Eh + L32L12)] ; 

aia = 5 sy (RE Ry T R; ois p= RR > 

BibT GG, = SE = SEE, DD) ; (D5) 
EA = Gs, EE = g5w) Ris] , 

d. = Buel Rh t Y3 En + ya Lo Ria + yr Los as) + A! Jg P ; 
9b s = Ga (aL Rh ur ye DR + Y3 TLR + Ua SPD RO) + de y^ DIRE , 
ad. m gs Rh t ys DS + ys LA Rio + 95 rn Ria) + gu SR ; 
Eu = Ge (aL Ra + ys DRE Tus LR JS P Lh Ry) + de y D ; 
Suus = a Us Ris + ys DF + yp Lo Ros y Dos Fs) + + ay! LL, 
ban =-2 lat R31 + D. L5,F, + ys Lo R32 + ya Lo RS) + E: y™ (GR i 


e vo" 
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The decay widths of the lightest new down-type quark b; are 


Mp, 
32T 
ID Ae — 22 LU feele Lie al 


W W 
+ 12:95, ne Date i 


T(bı > Wt) = ll + aiy + 0% — Any — 207 — Baier IT 


D, — Zb) = Z2: (1 c a5 + a} — 225 207 — 22207)" 


{(1 + a5 — 22 + (lee, + (Giada) | (D6) 


Z Z 
S 12295, rbr Ibinbr] ? 


m 
T(b, — h?b) = z4 + rfo a — 2x20 — =)? 


0 0 


h? h? 
+ ALIF, rbrIbrbin Ä ? 


where zr; = m¿/my, for index i = W,Z,b,h?. The heaviest new down-type quark bọ has 


six 


decay channels, the decay widths have the similar forms and can be deduced straight- 


forwardly. 
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